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ABSTRACT 

Angular momentum loss to circumbinary gas provides a possible mechanism for overcoming the 'last parsec' 
problem and allowing the most massive black hole binaries formed from galactic mergers to coalesce. Here, 
we show that if gas disks also catalyze the merger of the somewhat lower mass binaries detectable with the 
Laser Interferometer Space Antenna (LISA), then there may be a purely gravitational wave signature of the 
role of gas in the form of a small but finite eccentricity just prior to merger. Numerical simulations suggest that 
eccentricity, excited by the interaction between the binary and surrounding gas disk, is only partially damped 
during the final phase of gravitational radiation-driven inspiral. We estimate a typical eccentricity at one week 
prior to coalescence of e w 0.01. Higher terminal eccentricities, which can approach e = 0.1, are possible if the 
binary has an extreme mass ratio. The detection of even a small eccentricity prior to merger by LISA provides 
a possible discriminant between gas-driven inspirals and those effected by stellar processes. 

Subject headings: accretion, accretion disks — black hole physics — gravitational waves — galaxies: active 
— galaxies: nuclei 
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1. INTRODUCTION 

Supermassive black hole binaries formed from galactic 
mergers can coalesce if angular momentum loss - initially 
to sta rs and gas and sub sequently via g ravitational radia- 
tion ( Begelman, Blandford & Rees 1980) - is fast enough 
to conclude the process within at most a Hubble time, 
or more stringently before a subsequent merger creates a 
3 -bod y system that is unstable to a black hole ejection 
dValtonen et alJ[T994l) . Due to the formation of a loss cone, 
stars on their own are likely to fail to meet this require- 
ment for binary separation s of the order of r ~ 1 pc and bi- 
nary masses M > 10 8 Mrr, dMiiosaylievic & MerritfeOOltlYul 
120021 iBerczik. Merritt & Spurzemll2005l) . parameters which 
toda y are characteristic of massive (a ~ 200 km s" 1 ) galax- 
ies jGebhardtetafl |200(J |liniriie]fcMe^|2000). For 
this subset of galaxies, angular momentum loss to gas 
provides a plausible mechanism that can y i eld a prompt 
merger (iBeeelman. B landford & Rees 1980; Gould & Rixl 
l2000tlArmitage & Nataraianl2002HEscala et alJl2004l 120051 
At lower masses, stellar dynamic processes are in principle 
sufficient to effect mergers, but gas may still be involved if 
disks are present and remove angular momentum efficiently 
on a time scale significantly shorter than stellar interactions. 

Gas disks that catalyze binary mergers may generate ob- 
servable electromagnetic counterparts to the primary gravi- 
tational wave signature detectable with the Laser Interfer- 
ometer Space Antenna (LISA). Induced accretion or out- 
flows from small accretion disks surrounding the individual 
holes might yield precursors to the gravitational wave event 
(Armitage & Nataraian 2002), while the inflow of circumbi- 
nary gas left stranded as the binary contracts under gravita- 
tional radiation losses provides a robust prediction of an af- 
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terglow (Milosa vlievic & Phinnevl [2005). On a longer time 
scale, i mpulsive changes to the sp in of the hole following 
merger (Hughes & Blandford 2003) are likely to change the 
direction of any j ets launched from the immediate vicinity 
of the black hole (Merritt & Ekers 2002). Such counterparts 
are of interest, in part, because their detection in coincidence 
with LISA events would overcome LISA'S limited capability 
to spatially localize signals, and pin down the galaxy within 
which a merger had occurred. Unambiguous identification of 
merger sites significantly enhances the potential cosmological 
utility of merger events ( Holz & Hughes 2005; Kocsi s eit alJ 
120051) . and substantially improves the ability of LISA to trace 
the growth of black holes as a function of galactic environ- 
ment and morphology. 

In this paper, we study the evolution of binary eccentric- 
ity in the event that angular momentum loss to gas dom- 
inates the inspiral immediately prior to the onset of rapid 
shrinkage via gravitational radiation. In §2 we show that 
the gravitational interaction between the binary and a sur- 
rounding gas disk is likely to excite black hole binaries to 
eccentri cities e > 0.1, a res ult anticipat ed from bo th analytic 
studies (Artvmowicz 1992; Goldreich & Sari 2003), simula- 
tions of pre-main sequence binary stars tArtvmowi cz et alJ 
1991), and s imulatio ns of massive disk-embedded planets 
(Papaloizou. Nelson & Masset 2001). In §3 we argue that al- 
though gravitational radiation will damp these eccentricities 
prior to coalescence, the transition between disk-driven and 
gravitational wave-driven inspiral can occur at small enough 
radii that a small but significant eccentricity survives when the 
separation is only tens of gravitational radii. Detection of this 
eccentricity in the gravitational wave signal would provide 
strong circumstantial evidence for the role of gas in driving 
mergers, even in the absence of other electromagnetic coun- 
terparts. 

2. ECCENTRICITY DURING GAS-DRIVEN INSPIRALS 

If gas is the causative angular momentum loss process 
of supermassive black hole binary coalescences, the merger 
is expected to proceed through two ph ases. In the initial 
phase, the binary is either suit ounded by ( iGould & RixfeOOOt 
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lArmitage & Nataraianl2002l or completely embedded within 
JEscala et alJ l2004) circumbinary gas. If the gas mass is not 
too large, and the gas temperature (relative to the virial tem- 
perature at the radius of the binary) no t too high, gravitational 
torques (Goldreich & Tremaine 1980) fr om the binary act to 
clear a gap in the circumbinary disk dLin & Papaloizoui 1980). 
The inner truncation radius of this gap depends upon the bi- 
nary mass ratio, disk properties, and eccentricity. For small 
eccentricity truncation occurs at around 2a, where a is the bi- 
nary semi-major axis ( Artvmowicz & Lubowl [l994l) . The bi- 
nary, and the inner edge of the circumbinary disk, shrink as 
angular momentum is transfered from the orbital angular mo- 
mentum of the black holes into the gas. This phase is qualita- 
tively similar to the Type II migration (i.e. involving a gap) of 
massive planets, with the major difference being that a close 
black hole binary interacts with a relatively smaller mass of 
disk gas than a typical massive planet. The mismatch between 
the disk mass and the binary mass slows the rate of both gas 
inflow and orbital decay. 

The first phase ends when angular momentum loss to grav- 
itational radiation - which on its own drives decay of the 
semi-major axis a at a rate a oc a -3 - first exceeds that due 
to circumbinary disk torques. The binary then enters a sec- 
ond phase in which gravitational wave losses control the rate 
of inspiral. Initially, the circumbinary gas, although now dy- 
namically unimportant, flows inward fast enough to remain in 
contact with the tidal barrier created by the binary. Any initial 
eccentricity excited during the first phase will start to damp. 
Eventually, however, the decay of the binary separation be- 
comes too fast, and the binary detaches from the circumbi- 
nary gas, which through coalescence remains approximately 
frozen in the co nfigurati on it had at the moment of d etach- 
ment. (lArmitage & Nataraianl2002LlMilosavlievic & Phinnevl 
2005). Any low density gas remaining within the central cav- 
ity during this final stage is e xpected to have little influence 
on the binary orbital elements (Naravan 2000). 

2.1. Numerical simulations 

We have simulated the interaction between a binary and 
surroundin g circumbina ry gas disk using the Zeus code 
IStone & Normanlll992l) . We work in 2-dimensional cylin- 
drical co-ordinates (r, <p), with 'scaled gridding' in the radial 
direction (i.e. r,+i = (l + e)r,, with e a constant, so that com- 
putational zones have the same shape at all radii) and a stan- 
dard resolution of n r = 400, n^ = 400. The computational do- 
main extends from r m to r out , with r; n chosen to lie interior to 
the expected inner truncation radius of the circumbinary disk 
llArtvmowicz & Lubowll9 94). We exclude the region close to 
the black holes, and do not attempt to model the gas flow or 
torques arising from the smaller disks around each black hole. 
Outflow boundary conditions are imposed at both r m and r out . 

For a thin disk, the geometric thickness h jr is related to the 
local sound speed c s via h/r ~ c s /vk, where vk is the Keple- 
rian orbital velocity about the binary's center of mass. We 
consider self-similar disk models in which the equation of 
state is isothermal at any radius (P = pc 2 ), with an imposed 
radial scaling, 

c s (r)ocr- l/2 . (1) 

We adopt units in which G = 1, the sum of the masses of 
the black holes Mi +M2 = 1, and consider models in which 
h/r = 0.2 (corresponding to c s = 0.2 at r = 1), 0.1 and 0.05. 
We refer to these as the 'hot', 'warm' and 'cold' disk models 
respectively. 



At the relatively small radii we are most interested in (r < 
lOVj, where r s is the Schwarzschild radius), angular momen- 
tum transport within the circumbinary disk probably arises 
from tu rbulence driven by the Magnetorotational Instability 
(MRI) (Balbus & Hawlev 1991), which cannot be simulated 
in two dimensions. Instead we model angular momentum 
transport within the disk using a constant kinematic viscos- 
ity v, which operates on the azimuthal co mponent of the mo- 
mentum equation only (Papaloizou & Stanlev 1986). There 
is no reason to think that this, or any other simple viscos- 
ity prescription, reproduces accurately the multi-dimensional 
effects of angular m omentum transport from MHD turbu- 
lence (lO gilvid 120031) . We adopt it purely for numerical 
simpli city. The ma gnitude of v is set by reference to the 
IShakura & Sunvaevl (119731) viscosity prescription, 



We choose v such that the equivalent a = 0.1 at r = 2a, where 
a is the initial binary separation. Consistent with the choice 
of a constant viscosity, we set the disk surface density E(r) 
to be uniform in the initial state, so that the initial accretion 
rate M(r) oc vY, is constant with radius. The disk models with 
varying h/r are set up with the same initial surface density. 

The dynamics of the binary, including forces from the cir- 
cumbinary gas, are integrated using a simple second-order 
scheme that is adequate for the relatively short duration (t ~ 
10 2 orbits) of our simulations. Over these time scales, realis- 
tic surface densities of the circumbinary disk are small enough 
that they would lead to negligible changes in the binary semi- 
major axis a and eccentri city e. For example, a s teady-state 
Shakura-Sunyaev disk ( Shakura & Sunyaev 1973) around a 
black hole of mass 10 6 M Q , accreting at the Eddington limit 
with a = 0.1, has a surface density at 10 3 r s of the order of 
10 s gem" 2 . This means that the mass in the region of the disk 
that interacts with the binary is only a very small fraction of 
the binary mass, and evolution is slow. In order to see signif- 
icant evolution within of the order of 10 2 orbits, we therefore 
need to adopt a higher scaling for the surface density. How- 
ever, this does not affect the rate of eccentricity growth (or 
damping), measured as ade/da, because the perturbing accel- 
erations that alter a and e are linear in the disk mass (this is ev- 
ident explicitly in Gauss' equations, given in simplified form 
as equations (l)-(3) of Artymowicz et al. 1991). Choosing a 
higher surface density scaling therefore accelerates the time 
scale for changes in both a and e by the same factor. Large 
disk masses do lead to some wandering of the binary as the 
eccentricity grows - we suppress this and maintain the binary 
center of mass at r = throughout the runs. 

The initial conditions are set up as a uniform disk in Ke- 
plerian rotation (including the correction due to the pressure 
gradient) outside r = 2a, where a is the initial binary separa- 
tion. Inside 2a, we exponentially truncate the surface density. 
For the hot disk run, we evolve the disk for 40 binary orbits 
before turning on (smoothly over an additional 5 orbits) the 
backreaction from the disk in the binary integration. For the 
cooler disks, we start with the initial conditions for the hot 
disk and linearly reduce the sound speed over the first 25 or- 
bits of the calculation until it reaches its final, fixed value. 
These procedures are intended to minimize dynamical time 
scale transients created as the initially uniform disk adjusts to 
the non-axisymmetric forcing from the binary. The inner edge 
of the circumbinary disk, for example, is able to adjust its po- 
sition in response to the forcing on such a short time scale. 
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FIG. 1 . — Surface density at t = 40 orbits, just prior to starting active inte- 
gration of the binary, from the hottest disk with h/r = 0.2. The image is 8a 
on a side, where a is the initial binary separation. For this run the mass ratio 
q = Mi/Mi = 1 /3, and the initial eccentricity was e = 0.02. 



Considering longer time scales, however, the initial disk sur- 
face density profile - which corresponds to a steady-state so- 
lution around a single central mass - is not in equilibrium. 
Over a viscous time scale, gas whose inflow is impeded by 
tidal forces from the binary would accumulate near the inner 
edge of the circumbinary disk, increasing the surface density. 
Over the duration of our simulations, this effect is small but 
noticeable. It is not present for the hottest disk, which is able 
to overflow the gap toward the individual binary components 
most easily. 

2.2. Results 

Figure 1 shows a map of the surface density at t = 40 orbits 
for the hot disk, which surrounds a binary with mass ratio q = 
Mo/Mi = 1 /3, initial semi-major axis a = 1.02, and initial ec- 
centricity e = 0.02. The computational domain extended from 
r; n = 1.15 to r out = 20. Torques from the binary truncate the in- 
ner edge of the circumbinary disk at ar ound 2a, alth ough there 
is sig nificant overflow across the gap ( Artv mowicz & Lubowl 
1996) and through the inner boundary of the simulated disk. 
The cooler disk runs are qualitatively similar, although the 
thickness of the spiral arms and the extent of mass overflow 
both decrease as the disk becomes thinner. 

Figure 2 shows the evolution of a and e for binaries em- 
bedded within the cold, warm, and hot disks. Once the back- 
reaction forces from the disk are turned on after the initial 
equilibration period, the eccentricity grows and the orbit de- 
cays in all three cases. As expected, the rate of orbital decay 
is fastest for the hottest, most viscous disk model, while ec- 
centricity grows most rapidly for the coldest disk. This is con- 
sistent with the analytic expectation that eccentricity growth 
is driven by resonant interaction with gas at relatively large 
distance from the binary. For all three of our disk models, 
we find that the inner edge of the disk at the end of the equi- 
libration period remains close to r ~ 2a. The dominant ac- 
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FIG. 2. — Evolution of the binary eccentricity (upper panel) and semi-major 
axis (lower panel) for simulations of binaries embedded within disks of vary- 
ing geometric thickness. The eccentricity has been averaged over the orbital 
period to filter out small short time scale variations. The rate of eccentric- 
ity growth increases with decreasing disk thickness, while the rate of decay 
of the separation is reduced. Note that the normalization of the disk surface 
density has been chosen (arbitrarily) to yield measurable evolution over the 
~ 100 orbit time span of the simulations. 



tive resonance is then the 1:3 outer eccentric Lindblad reso- 
nance at r = 2.08a (Art vmowicz & Lubowi ri994). The rela- 
tive importance of this resonance is increased for cold disks, 
in which the inner evacuated cavity is better defined. For the 
cold disk model, although we are unable to determine the sat- 
uration value of the eccentricity, we find that e increases to 
more than e = 0.1 before the binary has shrunk by more than 
a few percent in semi-major axis. 

Figure 3 shows the evolution of binary orbital elements in 
the a-e plane, which for the reasons noted above is indepen- 
dent of the arbitrary choice of disk surface density. For the 
cold disk we find that the eccentricity grows at a relatively 
constant rate, 

de 

a— ~ -2.05, (3) 
da 

up to e w 0. 1, with the growth rate modestly decreasing as the 
eccentricity increases further. Any significant orbital migra- 
tion within such a disk would therefore be expected to yield 
significant eccentricities. Growth rates are substantially lower 
in the warm (ade/da ~ -0.25) and cold (ode/da as -0.1) 
disks. This trend is qualitatively consistent with how one 
might expect da / dt and de / dt to scale with binary parameters. 
Defining the ratio of the local disk mass to the secondary mass 
via qa = ira 2 'E/M2, we expect the disk interaction to shrink the 
binary orbit at a rate, 



da 

— oc —a 
dt \r 



q d Qa. 



(4) 



This is just the viscous inflow velocity of gas in the disk, re- 
duced by the ratio of the local disk mass to the secondary 
mass. We note that this expression does not include changes 
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FIG. 3. — Binary evolution in the a-e plane for disks of varying geometric 
thickness. The highest initial rates of eccentricity growth (ade/da ~ —2), 
shown as the dashed line) are obtained for the coolest, thinnest disk simulated. 
The arrow shows the estimated growth rate for the even colder disk model 
simulated by Artymowicz et al. (1991). 



to the equilibrium disk structure that occur over a viscous time 
scale dSver & ClarkdH995l) . Viscous time scale equilibrium, 
which is not achieved in our simulations, is discussed more 
thoroughly in §3. The rate of change of ecce ntricity is ex- 
pected to scale as (Gold reich & Tremainell9 80). 

de n 

— oc q qMe. (5) 

Comparing these expressions we expect eccentricity growth 
to be most rapid for nearly equal mass ratio binaries in geo- 
metrically thin and / or relatively inviscid disks. 

2.3. Comparison with previous work 

Our simulation s are similar in spirit to those of 
lArtvmowicz et alJ (11991ft and iGunther. Schafer & Kiev I 
(2004 ). Using smooth p article hydrodynamics simulations, 
Arty mowicz et alJ {l99 1 ) derived an eccentricity growth rate 
ade/da « -4.4 for a binary with a mass ratio of q = 3/7 
and e = 0.1, embedded within a disk in which h/r ss 0.03 
and a ~ 0.1. This disk is modestly colder than our 'cold' 
case. Their growth rate, shown as an arrow in Figure 3, is 
approximately a factor of two greater than obtained in our 
'cold' simulation. Our smaller growth rate is consistent with 
a trend toward faster growth in geometrically thinner disks. 

Motivated by the surprisingly large eccentricities of many 
massive extrasolar planets, more recent numerical work on 
binary-disk interactions has focused on extreme mass ratio 
(q <C 1) systems. In a disk model intend ed to represent a typi- 
cal protoplanetary disk, Papaloizou. Nelson & Masset (2001) 
obtained eccentricity growth for mass ratios q > 0.02, while 
damping occurred for more extreme mass ratios. Eccen- 
tricities as large as e w 0.25 were obtaine d in these simu- 
lations. Su bsequent analytic work (Oailvie & Lubowl l2003t 
lGoldreich&Saril2 003) has identified circumstances in which 



the critical mass ratio for eccentricity growth might be re- 
duced still further - perhaps to values of q ~ 10~ 3 - but this 
has yet to be confirmed numerically. In the black hole bi- 
nary problem, these subtleties are probably moot. Extreme 
mass ratios are unlikely to occur as a consequence of galac- 
tic mergers llYul 120021) . because the dynamical friction time 
scale for the smaller black hole to reach the galactic center 
would exceed the Hubble time. Extreme mass ratio bina- 
ries are however conceiv able if intermediate mass black holes 
( iMiller & Colberi2004l) form within the nucleus and are sub- 
sequently captured by the gas disk. 

2.4. Initial eccentricity 

Our simulations commence with finite eccentricity, which 
is then unstable to subsequent growth. The physical ori- 
gin of such a seed eccentricity is speculative, although there 
are enough possibilities that we think it unlikely that a bi- 
nary could evade the instability by remaining perfectly cir- 
cular. First, the endpoint of the prior stage of stellar dy- 
namical hardening might leave the binary with non-zero ec- 
centric ity, although the results of Berczik. Merritt & Spurzem 
( 120051) suggest that any such eccentricity would be small. 
Second, the binary could acquire eccentricity while embed- 
ded within the disk, for example via a Kozai resonance with 
a misaligne d gas or stellar dis k at larger radius in the nucleus 
( iBlaes. Lee & Socrate s 2002). Finally, non-axisy mmetry in 
the disk itself - due to self-gravity at large radii (Goodman 

f 03|). eccentr ic instabilities in individual circumstellar disks 
ubowl [l99l|). or in stabilities at the edge of the circumbi- 
nary disk (Goldreich & Sari 2003) - appears to be almost un- 
avoidable. Such non-axisym metry in the disk is readily ex- 
changed into the binary IPapaloizo u. Nelson & Massell2001t 
Papaloizou 2002), and, if it is large enough, could serve as the 
initial seed for subsequent growth. 

3. ECCENTRICITY OF BINARIES CLOSE TO COALESCENCE 

Any eccentricity attained by the binary during the phase of 
gas-driven orbital decay will be damped during the final stage 
of gravitational wave-driven inspiral. The extent of the damp- 
ing depends upon the separation of the binary (in gravitational 
radii) at the time when decoupling from the gas occurs. To 
gain a rough idea of the magnitude of the damping, we ini- 
tially consider a toy model in which the binary decouples in- 
stantaneously from the disk at a radius a cr ; t , which for now 
we regard as a free parameter. For a < a cr ; t , the binary orbit 
evolves solely under the influence of gravitational radiation 
losses. For a bi nary with energy E an d angular momentum /, 
at such a stage (Padmanabhan 2001a), 

dE _ 32G 4 M 2 l Ml(M l +M 2 ) 1 
dt ~ 5c 5 a 5 (l-e 2 )i 

73 37 

L 24 96 ' 
dJ _ 32G2 A%M%s/Mi+M 2 1 
dt ~ 5^1 X (l-e 2 ) 2 

x[l+ 7 -e 2 ]. (7) 

These equations are readily integrated to yield the resid- 
ual eccentricity at a fiducial time (typically either one week 
or one year) prior to coalescence, as a function of a cr i t , 
q, and the assumed initial eccentricity generated during the 
earlier gas-dominated phase. The formulae apply in the 
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FIG. 4. — Binary eccentricity at a time to prior to coalescence in a toy 
model in which the binary is assumed to attain an eccentricity of e = 0.25 
during the initial phase of gas-driven inspiral at a > a crit . At a ait we assume 
that gravitational radiation starts to dominate over angular momentum loss to 
the gas, so that for a < a a - It inspiral and attendant damping of e is governed 
solely by gravitational radiation losses. The solid lines show results for binary 
mass ratios q = Mi/Mi of 1, 0.1, 0.01 (the upper solid curve is for q = 1), 
and to = 1 week. The dashed curves are for to = 1 year. The total system 
mass is 10 6 Mq in all cases. 



weak field limit, which is appropriate provided that we do 
not consider times extremely close to the instant of coales- 
cence (note that at one week before coalescence, the sep- 
aration is still a w 4-QGM/c 2 ). In the strong field regime, 
e tends to increase in anticipation of the final 'plunge' 
JCufler. Kennefick & Poissonl 19941) . 

Figure 4 shows the predicted binary eccentricity at one 
week and one year before coalescence for different q, as 
a function of the decoupling radius a cl -i t . An initial eccen- 
tricity e = 0.25 was assumed, whic h seems r easonable given 
our simulations and previous work lArtymowicz et alj| 199U 
lArtvmowiczll 1 9921 iPapaloizou. Nel son & Masset EOOll) . Bi- 
naries observed through coalescence with LISA are predicted 
to have a mass ratio distribution that peaks around q ~ 0.1. 
If these binaries in fact have a gas-driven inspiral phase, de- 
coupling within 10 3 gravitational radii typically yields a final 
eccentricity of e ~ 0.01. At fixed a cl -i t , binaries with more ex- 
treme mass ratios yield smaller final values of the eccentricity, 
reflecting the fact that at a given time prior to coalescence the 
separation of q <C 1 systems is smaller than for the q = 1 case. 

To go further, we need to estimate how a cr j t depends upon 
the binary and disk parameters. We begin by noting that, in 
the absence of a binary, gas in an accretion disk flows inward 
at a rate given by the viscous inflow speed, 



3 fh 



dv\s.c. — _ I 



av K . 



(8) 



The separation of a circular binary shrinking under the influ- 
ence of gravitational radiation decays at a rate, 



flGW : 



64G 3 MiM 2 (Mi +M 2 ) 



(9) 



Equating these expressions determines a detachment ra- 
dius, defined via a'visc(adetach) = a'GwOdetach). The de- 
tachment radius is the smallest radius for which the in- 
flowing gas can 'keep up' with the shrinkage of the bi- 
nary. Interior to this radius, the binary orbit continues 
to contract within a growing hole while the circumbinary 
gas remains essentially frozen jArmitage & Nataraianl l2002t 
iMilosavlievic & Phinnevl 120051) . The detachment radius is 
typica lly quite small. For example. IMilosavlievic & Phinnevl 
(2005) estimate that the shrinking binary leaves a hole of size 
i=s 120GM/c 2 , which lies well within the innermost, radiation 
pressure dominated disk. 

For dynamical purposes we are more interested in the de- 
coupling radius a cnt - within which the rate of angular mo- 
mentum loss to gravitational waves exceeds that due to inter- 
action with the circumbinary disk. To calculate this we need 
to set q'gw equal to the disk-driven inspiral velocity ai nsp irai- 
The inspiral velocity equals the viscous inflow speed in the 
limit where the secondary black hole is of low enough mass 
that it can be treated as a test particle within the disk. If the un- 
perturbed disk surface density at the orbital radius r of the sec- 
ondary is So, test particle behavior will occur for M 2 < 7rr 2 Eo- 
For more massive secondaries, the large reservoir of binary 
angular momentum reduces the inspiral velocity to a fraction 
of the viscous inflow speed. The binary then acts as a slowly 
moving dam which impedes the free inflow of gas. 

For black hole binaries the massive secondary limit is the 
relevant one for small enough separations. To estimate the 
inspiral veloci t y in t his regime, we apply the analysis of 
iSver & Clarkd i 19951) . For a di sk in which the unpe rturbed 
surface density at radius r is £o,|S}fer & Clarke ( 1995) define 
a measure of disk dominance, 



B = 



47rr 2 S 



(10) 



where M 2 is the secondary mass. Then if the unperturbed disk 
solution can be locally written as, 



«#r 7 , 



(11) 



(12) 



with j3 > 0, the inspiral velocity is, 

^'inspiral — ^'visc ^ 1) 
flinspi ra l = flviscB /3/(1+/3) (B<1). 

The correction factor in the last equation accounts for the 
slowdown in orbital decay when the secondary is massive 
compared to the local disk. The slowdown is usually less than 
the full factor of B, because gas accumulates just outside the 
slowly moving tidal barrier created by the binary. This pile-up 
increases the torque rela tive to tha t obtained from an unper- 
turbed disk solution (Pringle 1991). 

The above equations allow us to calculate the decoupling 
radius for a binary given an estimate for a and a model for the 
accretion disk structure. Here, we adopt the simplest Shakura- 
Sunyaev (1973) model. We normalize the accretion rate to the 
Eddington limiting value, 



M 



Mi 



Edd 



M Edd = 2.2 x 10- 



M 



M yr" 



(13) 
(14) 



and measure the mass m = M /Mq in Solar units. Then in the 
inner, radiation pressure-dominated disk, away from the inner 
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boundary ( Padmanabhan 2001b), 



— =7A6m, 
r. 



E = 0.42a" 1 m" 1 (r/r. s ) 3/2 gem 



(15) 
(16) 



whereas at larger radii, where gas pressure dominates and the 
opacity is due to electron scattering, 



-l.exlO-V'/'V'^m-'/'V/r,) 21 



/20 



(17) 



S = 9.7 x 10V 4/5 ;7i 3/ V /5 (rA. s r jp gcm" z . (18) 

Here r s = 2GM/c 2 is the Schwarzschild radius, and the transi- 
tion between the two regimes occurs at, 



ii = 360a 2 / 21 m 16 / 21 m 2 / 21 . 
r. 



(19) 



These relations allow us to calculate a' v isc in both the radiation 
pressure and gas pressure dominated regions of the disk. In 
the gas pressure dominated middle disk we can also directly 
estimate a'mspirai- In this region (3 in equation (lilt is 0.6, so 
the slowdown in the inspiral velocity when the secondary is 
massive scales as B° . For representative disk and binary 
parameters (M l =9x 10 5 M , M 2 = 10 5 M Q , a = 10" 2 , m = 1) 
we find that the disk drives inspiral at a rate reduced from the 
viscous rate by about a factor of 50. We cannot use the same 
formula to calculate the slowdown in the radiation pressure 
dominated inner disk, since in this region j3 < 0. For an esti- 
mate, we simply continue to use the slowdown factor derived 
for the middle disk in this region. We note that the struc- 
ture of accretion disks in the innermost radiation-dominated 
regim e remains su bject to large uncertainties (iTurneJ 120041: 
iTurner et al.l2005l) . so any statements about this region are at 
best semi-quantitative. 

The solid curves in Figure 5 show our best estimates for 
the estimated eccentricity at one week and one year prior to 
coalescence. We take a = 0.01, m = 1, and consider binaries 
of varying q but fixed total mass Mi +M 2 = 10 6 M Q . Calcu- 
lations suggest that LISA should detect a significant number 
of mergers in which the black holes have masses in this range 
( Sesa na et al.ll200 5). whereas events involving more massive 
holes (M ~ 10 8 Mq) will be rare. As previously, we assume 
that e = 0.25 at the epoch when gravitational radiation be- 
comes the dominant angular momentum loss mechanism and 
the binary decouples from the gas. With these assumptions, 
we find that for roughly equal mass binaries e ~ 7 x 10~ 3 at 
one week before coalescence, while a year earlier e ~ 0.02. 
Higher terminal eccentricities are predicted for more extreme 
mass ratios, first because gravitational radiation is relatively 
less efficient, and second because the slowdown in disk driven 
inflow at small radii is smaller for a less massive secondary 
(eauationll2>. The eccentricity at one year before coalescence 
could exceed e = 0.1 for q < 0.02. We also note that there is a 
substantial separation between the decoupling radius - when 
the gas ceases to be dynamically important for the evolution 
of the binary - and the detachment radius when the gas and 
the binary become physically separated. 

With these parameters the derived decoupling radius a cr i t 
moves from modestly interior to r t (for q = 1), to well inside 
the radiation pressure dominated disk (for q <C 1). The un- 
certainties in calculating a'inspirai in this region therefore have 
a direct impact on the derived values of a cr i t , and on the fi- 
nal estimated eccentricity. Cognizant of this, we also plot in 



o.i 



0.01 



0.001 




FIG. 5. — The predicted binary eccentricity at 1 year (upper solid curve) 
and 1 week (lower solid curve) prior to coalescence as a function of mass 
ratio q. We assume that m = 1, a = 0.01, set the initial (gas-driven) eccen- 
tricity to e = 0.25, and the total system mass M t +M 2 = 10 6 M Q . With these 
parameters, the transition between disk-driven and gravitational wave-driven 
orbital decay occurs close to the transition radius in the disk between radia- 
tion and gas pressure domination (for q = 1), and moves to smaller radii for 
q <^.\. Since the migration rate in the radiation pressure dominated disk is 
uncertain, we also plot a lower limit to the final eccentricity (dashed curves), 
obtained by assuming that gravitational radiation always dominates in the ra- 
diation pressure dominated disk. The dashed curves are evaluated at 1 year 
and 1 week prior to merger. 



Figure 5 a more robust lower limit to the final eccentricity, de- 
rived under the simple assumption that the disk dominates the 
loss of angular momentum up to the point where it becomes 
radiation pressure dominated. We assume an unperturbed disk 
structure when determining this radius (equationll9>. so there 
is no dependence on q. For q > 0.2, the resulting lower limit 
is approximately e ~ 0.01 at one year before coalescence, and 
e~3x 10~ 3 at one week before coalescence. Since the de- 
coupling radius for this calculation is fixed, the trend is for 
smaller final eccentricity for lower q (e.g. Figure 4). 

4. DISCUSSION AND SUMMARY 

In this paper, we have argued that observation of the gravi- 
tational waves from merging black hole binaries could pro- 
vide circumstantial evidence of the uncertain astrophysical 
processes that (probably) allow supermassive black hole bina- 
ries to coalesce. Specifically, we suggest that a small but non- 
zero eccentricity would signal that angular momentum loss to 
gas was the dominant decay process immediately preceding 
the final gravitational wave-driven inspiral. If, alternatively, 
interaction with stars precedes the final inspiral, the final ec- 
centricity is likely to be very small. Most studies suggest 
that orbital decay via st ellar dynamical processes is accompa- 
nied by circularizatio n ( Polnarev & Rees 1994; Makin oIl997t 
lOuinlan & Hernauistll 19971) . although I Aarsethl J2003I) has re- 
ported eccentricity growth in some recent N-body simula- 
tions. 

In the case of gas-driven mergers, eccentricity 
growth appears to be a robust prediction of simulations 
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jArtvmowicz et all 1199 It iPaoaloizou. Nelson & M asset 
12001 " provided that (a) the mass ratio is not very extreme 
(q ~ 0.1 appears to be safely in the 'growth' regime), and 
(b) the surrounding gas is not so hot and dense as to com- 
pletely swamp the binary and thereby fill in th e gap whose 
existe nce is necessary for eccentricity growth ( Escala et al. 
2004). One caveat is that existing simulations do not treat 
the angular momentum transport within the disk from first 
principles, and it is possible that this omission might alter the 
conditions needed for eccentricity growth. Direct simulations 
of binaries embedded within MHD turbulent disks are 
computationally possible (altho ugh still difficult) only in the 
gas pressure dominated regime jNelson & Papaloizou 2003; 
Winters, Balbus & Hawlev 2003), and would be valuable. 
There is even more uncertainty as to how the binary interacts 
with a radiation pressure dominated disk. 

The presence of significant eccentricity for black hole bi- 
naries close to merger has additional consequences that are 
potent ially observable. As noted by Artvmowicz & Lubow 
( 1996) both in general, and for the specific case of the quasar 
OJ 287, accretion from the circumbinary disk on to the in- 
dividual disks around the holes is strongly time-variable if 
the binary is eccentric. The time scales can be interesting 
even when the binary is still far from coalescence. For exam- 
ple, a circumbinary disk surrounding a binary with M\ +M2 = 
10 6 M Q and q = 1 /4 would remain gas pressure dominated un- 
til the separation was less than w 860r s , at which point we 
would predict the disk torque to dominate and the eccentricity 
to be significant. The inspiral time for these parameters is of 
the order of 10 5 yr, but the orbital time scale is still small - 
approximately 1 month. This is also the time scale on which 
the individual disks might be expected to vary. If gas is re- 
sponsible for merging a large fraction of supermassive black 



hole binaries, such systems, because they are that much fur- 
ther away from merger, must be much more numerous than 
galaxies that will host LISA events. Although still rare, they 
might be identified via their characteristic X-ray variability. 

The eccentricity is one of the intrinsic parameters of merg- 
ing black holes that can be measured via detection of the 
gravitational wave emission, potentially to high accuracy. As 
an illustration, for a stellar mass black hole merging with 
a 10 6 Mq supermassive black hole from an orbit with e = 
0.1, LISA may achieve an accurac y of Ae ~ 10" 4 for a n 
event with signal-to-noise ratio of 30 (Barack & Cutler 2004). 
Mille r^tafl d20 05 ) discuss how LISA measurements of the 
eccentricity of compact objects merging with supermassive 
black holes can be used to distinguish between different cap- 
ture mechanisms. Similar measurements for mergers of two 
supermassive holes may allow the role of gas in the evolution 
of binary black holes to be studied as a function of the masses 
of the holes involved. 
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